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^ ■ Abstract: In the framework of semi-hard QCD approach, we present a consistent 

analysis of D* meson production at HERA energies. The consideration is based on 

^ , universal unintegrated gluon densities, which have BFKL behavior in the small x 

region. Predictions of the CCFM evolution equation for D* production are obtained 

CN I and show a good description of D* data at HERA. 

^ : 1 Introduction 
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The standard parton model is based on the DGLAP [ |T], ^ ^ Q evolution equations, which 
QJ ', re-sums contributions from [asln(yU^/A^)]. It represents an one-dimensional phase space ap- 

"^ ' proximation for the parton motion, also known as the collinear approximation, which gives 
the correct behavior of the structure functions at not too small values of x. When x becomes 

^ ' smaller, also contributions from [asln(/i^/A^) ln(l/x)] and [Q;sln(l/x)] need to be considered. In 
the so called kt factorization or semi-hard approach (SHA) [|^, ^, 0, P], the transverse momenta 
of the partons in the evolution from large x at the proton vertex towards small x at the hard 
interaction vertex are taken into account. This evolution in x has been formulated in terms 
of the BFKL [ I, |T|, [ill evolution equation. The CCFM [ |l|, |l|, 0, |15| evolution equation 



S 



includes coherence effects via angular ordering and it reproduces the BFKL (DGLAP) evolution 
equation in the small (large) x limits, respectively. 

The resummation [i, i, 0, H of the terms [asln(/i^/A^)], [Q;sln(/i^/A^) ln(l/x)] and [Q;5ln(l/x)] 
in SHA results in the so called unintegrated gluon distribution jF(x, k^, Ql), which determines 
the probability to find a gluon carrying the longitudinal momentum fraction x and transverse 
momentum kt. The factorization scale Qo (such that as{Ql) < 1) indicates the non perturba- 
tive input distribution. They obey the BFKL [ |, 0, [TTl or CCFM [ 0, |13|, |T|, [l|] equation 



and reduce to the conventional parton densities -F(x, /i^) once the kf dependence is integrated 
out: 

r ^(x, kl Ql) dkl = X F{x, fi', Ql). (1) 

Jo 

However, in CCFM the unintegrated parton distribution A{x, k^, Ql, (f) (instead of jF(x, k^, Ql)) 

depends also on the maximum angle for any emission corresponding to q (coming from angular 

ordering) . In the small x limit they reduce to JF [ [U 



To calculate the cross section of a physical process, the unintegrated functions J^ or A 
have to be convoluted with off-mass shell matrix elements [ 0, H corresponding to the relevant 
partonic subprocesses. In off-mass shell matrix element the virtual gluon polarization tensor is 
taken in form of SHA prescription [ ^: 

Llf^ = 4^ = P>V/|A;,p = k^k-^/lhl^ (2) 

The specific properties of semi-hard theory may manifest in several ways. With respect to 
inclusive production properties, one obtains an additional contribution to the cross sections 
due to the integration over the kf region above /i^ and the broadening of the pt spectra due 
to extra transverse momentum of the interacting gluons [ |^, |^, [16|, [l^ . It is important that 
the gluons are not on-mass shell but are characterized by virtual masses proportional to their 
transverse momentum. This also assumes a modification of the polarization density matrix. A 
striking consequence of this fact on the J/ip spin alignment has been demonstrated in [ |l^ . 

In this paper we present predictions for the production of D* mesons in photo-production at 
HERA using the SHA approach. We use an unintegrated gluon density coming from a solution 
of the CCFM equation (see [ ^). We also show predictions based on a parton level Monte 
Carlo integration using a BFKL-like parameterization of the unintegrated gluon density, and 
we compare both predictions. 



2 Unintegrated gluon distribution and CCFM evolution 

The parton evolution at small values of x is believed to be best described by the CCFM evolution 
equation [ |12], |13|, |l^, [l^, which for a; — > is equivalent to the BFKL evolution equation [ 



^ |Ty, |Tl| and for large x reproduces the standard DGLAP equations. The CCFM evolution 
equation takes coherence effects of the radiated gluons into account via angular ordering. In [ 
^] it is shown that a very good description of the inclusive structure function F2{x,Q'^) and 
the production of forward jets in DIS, which are believed to be a prominent signature of small 
X parton dynamics, can be obtained from the CCFM evolution equation. The main important 
point there was the treatment of the non-Sudakov form factor, which suppresses radiation at 
small values of x. In Fig. |l] we show the gluon density obtained from this solution of the CCFM 



equation [ |T9| as a function of x for different values of k^ at g^ = 10 GeV^. In Fig. |^ we show 



the gluon density as a function of k"^ for different values of x at g^ = 10 GeV^. 

For comparison, we also use the results of a BFKL-like parameterization of the unintegrated 
gluon distribution J-'{x, k^, fi"^), according to the prescription in [ ^]. The proposed method lies 
upon a straightforward perturbative solution of the BFKL equation where the coUinear gluon 
density xG{x,fi'^) from the standard GRV set [ ^ is used as the boundary condition in the 
integral form (0). Technically, the unintegrated gluon density is calculated as a convolution of 
coUinear gluon density with universal weight factors [ EOl : 



^(x, kl /i^) = /' g{r^, kl ^2) ^ G{-, /i^) dv, (3) 

Jx 7] rj 



as 



Giv^kl^^') = ^ Jo(2v'«.ln(l^)ln(/xVfcf)), k', < /x^ (4) 
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Figure 1: The gluon density A{x^k1,(f) (solid line) obtained from a solution of the CCFM 
equation [\T^] as a function of x for different values of k^ (at (f = 10 CeV^). Also shown for 
comparison is the unintegrated gluon density function jF(a;, fc^,/i^) according to the parameteri- 
zation of Bliimlein [^^ (dashed line) at fj,'^ = 10 CeV^ and of Kwiecinski, Martin, Stasto [^1\] 
(dotted line). 
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(5) 



where Jo and Iq stand for Bessel functions (of real and imaginary arguments, respectively), 
and ois = Sa^/vr. The latter parameter is connected with the Pomeron trajectory intercept: 
A = Q;s41n2 in the LO and A = as41n2 — A^a^ in the NLO approximations, respectively, where 

In the following we use A 



A^ is a number 23 



0.35. 



The presence of the two different parameters, iJ? and A;^, in eq.(^ for unintegrated gluon 
distribution JF(x, k^, /i^) refers to the fact that the evolution of parton densities is done in two 
steps. First the DGLAP scheme is applied to evolve the structure function from Ql to /x^ within 
the coUinear approximation. After that eqs.(|^)-(|^) are used to develop the parton transverse 
momenta k^. This is in contrast to the CCFM evolution, where the evolution of "longitudinal" 
and "transverse" components occurs simultaneously. 

From Figs. |1] and we see that the BFKL approach gives much harder kf spectrum than 



the CCFM approach. However, it has been argued extensively in the literature [ 21, ^ 25 
that in BFKL a so-called "consistency constraint" should be applied, so simulate at least a 
part of the large next-to-leading corrections. For comparison we also show in Figs. |l| and |] 
the unintegrated gluon distribution from Kwiecinski, Martin, Stasto [ pT|]p|. This shows that 
the shape of the distributions from BFKL including the "consistency constraint" is similar to 



^A. Stasto kindly provided us (H.J.) with the program. 
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Figure 2: The gluon density A{x,k'^,(f) (solid line) obtained from a solution of the CCFM 
equation [\T^] as a function of kl for different values of x (at (f = 10 CeV^). Also shown for 
comparison is the unintegrated gluon density function jF(x, A;^,/i^) according to the parameteri- 
zation of Bliimlein [^^ (dashed line) at fj,'^ = 10 CeV^ and of Kwiecinski, Martin, Stasto [^1\] 
(dotted line). 



the one obtained from CCFM and that the gluon density is stronger suppressed at large kt as 
compared to the approach in [ |2^ . Unfortunately we could not use the the gluon distribution 



from [|21[], because it started only at kf > 1. 



3 Predictions for D* meson production at HERA 



We have used the hadron level Monte Carlo program Cascade described in [ |T9[ to predict the 
cross section for D* photo-production at HERA energies. The unintegrated gluon distribution 
was obtained from the solution of the CCFM equation described in [ |T9[. The scale in as was 
set to k'f in the parton evolution and we used A^qcd — 0-2 GeV. For the hard scattering the 
off-shell matrix element for heavy quarks (including the heavy quark mass with nic = 1.5 GeV) 
as described in [ |^ are used together with the one-loop expression for as with k^ (of the gluon 
entering the hard scattering) as the scale. The complete initial state cascade is simulated via a 

I^. The hadronization was performed by the Lund string 
. The Peterson function with e = 0.06 was used for the 



backward evolution as described in 
fragmentation JETSET [ |6|, |3, | 



charm quark fragmentation. 

In Fig. ^ we show the prediction of D* production as a function of the transverse momentum 
pf* using the Cascade Monte Carlo described above and compare it with the measurement 
of the ZEUS collaboration [ ^ . We observe a rather good description of the pt spectrum. 
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Figure 3: The differential cross section da/dpf for Q^ < 1 GeV^. The solid line shows the 
prediction from the Cascade Monte Carlo and the dashed line is the calculation of [^9^. The 



data point are from [ p^ . 



In Fig. ^ we show the D* cross section as a function of the pseudo-rapidity r]^ for different 
regions in pt. Also here we observe a good description of the experimental data points. Here 
the main important point is the cross section at values oi rj^ > 0.5. 

For comparison we used the results of a calculation of [ ^ with a parameterization for the 
unintegrated gluon distribution jF(a;, fc^,/i^), according to the prescription of [ ^. Here the 
scale m^ is used in a^ The results are also shown in Figs. |^ and ^. In general both prediction 
agree rather nicely. The differences observed are entirely due to the different behavior of the 
unintegrated gluon distribution as a function of x and k"^. In addition in the CCFM approach 
angular ordering and the maximum angle allowed for any emission plays a important role. The 
results presented here are similar to the ones obtained from a full NLO calculation. In the SHA 
approach the gluon entering the hard interaction is off-mass shell, which is a similar situation 
as in a NLO calculation where the propagators in the 3 parton final states are fully considered. 
This is in contrast to the LO DGLAP (coUinear) case, where the gluon is always treated on- 
mass shell. However if, in LO DGLAP, heavy flavor excitation is included (via resolved photon 
processes) then again a similar situation occurs. It was found in [ ^ that including heavy 
flavor excitation in LO Monte Carlo programs, lead to better description of the data. The 
semi-hard approach presented here shows, that a good description of also the photo-production 
of heavy flavor can be achieved in a theoretically consistent way, without including artificially 
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Figure 4: T/ie differential cross section da/drj^' for Q^ <1 GeV^ for different regions of pf . 
The solid line show the prediction from the Cascade Monte Carlo and the dashed line is the 
calculation of f^J. The data point are from [13^]. 



large intrinsic transverse momenta, or heavy flavor excitation. 



4 Conclusions 



We have shown that using a unintegrated gluon distribution obtained from a solution of the 
CCFM evolution equation, that describes the structure function F2(x, Q^) and the production 
of forward jets in DIS at HERA, we can also describe the cross sections of inclusive D*^ meson 
production measured at HERA. Within the semi-hard approach the measured cross section as 
a function of pt and rjo* can be nicely described. The results are similar to NLO calculations. 
The shape of the gluon kt distribution is driven by the BFKL or CCFM evolution equations. 
This shows that there is no place to include any artificially large intrinsic transverse momentum 
distribution of parton inside the proton. 

It is also interesting to note, that within the semi-hard approach, heavy flavor excitation in 
the photon is consistently included, by the fact that a gluon radiated close to the quark box 
can have a transverse momentum larger that of the quarks. 
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